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ABSTRACT 

Weak gravitational lensing provides a means of testing the long-range properties of gravity. Current 
measurements are consistent with standard Newtonian gravity and inconsistent with substantial mod- 
ifications on Mpc scales. The data allows long range gravity to deviate from a 1/r potential only on 
scales where standard cosmology would use normal gravity but be dominated by dark matter. Thus, 
abnormal gravity theories must introduce two fine-tuning scales - an inner scale to explain flat rotation 
curves and an outer scale to force a return to Newtonian gravity on large scales - and these scales must 
coincidently match the scales produced by the dark matter theory after evolving the universe for 10 
billion years starting from initial conditions which arc exquisitely determined from the cosmic microwave 
background. 

Subject headings: cosmology: theory - gravitational lensing 



1. INTRODUCTION 

Weak lensing of background galaxies by foreground 
large-scale structure offers an opportunity to directly probe 
the mass distribution on large scales over a wide ra nge o f 
redshifts. As first poi nted out by Blandford et al. (1991) 
and Miralda-Escude 



1991 



these effects are of order a 
few percent in adiabatic cold dark matter models making 
their observation challenging but feasible. Early predic- 
tions for the power spe ctrum of the shear and convergence 
were made by Kaiser ( |1992| ) on the b asis of linear pertur- 
bation theory. Jain & Seljak ( [1997 ) estimated the effect 
of non-linearities in the den sity t hrough analytic fitting 
formulae (Peacock & Dodds 1996| ) and showed they sub- 
stantially increase the power in the convergence below the 
degree scale. Because weak lensing can measure the mat- 
ter power spectrum without many of the problems of ap- 
proaches based on the distributions of galaxies or clusters 
(e.g. bias), it may ultimately provide as clean a cosmolog- 
ical probe as the microwave background. Recently, several 
observational groups have report ed co nvincing evide nce of 
the effect (v an W aerbeke et al. |2000 ; Bacon et al. 2000 
Kaiser et al. |2000t Wittman et al. pOOq ; Maoli et al. pOOl 
Rhodes et al. |2001|; van Waerbeke et al. |200l|) 



All these theoretical and observational studies are pri- 
marily motivated by standard theories of gravity and cos- 
mology. Despite the tremendous overall success of these 
theories, there has been a recent resurgence of inter- 
est in non-standard theories of gravity, largely motivated 
by the possibility that the standard paradigm has diffi- 
culty matching the dyn amical st ructu re of galaxies (e.g. 
Flore s fc P rimack |1994 ; Moore |1994| ; Navarro & Stein- 
metz 2000| ). Most of these proposed modifications aim to 
make gravity a longer-ranged force on scales comparable 
to the sizes of galaxies in order to explain the flat rota- 
tion curves of galaxies on scales larger than the app arent 
distribut ion o f matter (e. g. Sellwood fc Kosows ky 2000, 
Sanders |l99gl |1999| , |2000| ; McGaug h |1999| , ^000| , but see 
Scott et al. 2001 and Aguirre et al. 2001 for an op posite 



perspec tive). As has been noted befor e (e.g . Krishcr 1988 
Walker J994; Bekenstein & Sanders |1994; Z hytni kovfe 



& Bernardeau 2001; Mortlock & Turner 2001) any longer 
ranged gravitational force, if it also affects photons, should 
have implications for gravitational lensing. In particular 
it should profoundly affect the strength of weak lensing 
shears on large scales. Many of the above authors, how- 
ever, consider gravitational lensing by isolated objects. To 
understand the lensing effects of modifying gravity on large 
scales it is necessary to use the weak lensing formalism, 
summing over the contributions from all density pertur- 
bations. 

2. THE MODEL 
We base our models on the discussion by Zhytnikov & 



Nester (1994) of modified gravity theories wi thin t he con- 
text of linearized relativity (see also Edery 1999[ ). This 
framework provides a relativistic gravity model which au- 
tomatically obeys the equivalence principle and within 
which definite calculations can be made, while at the same 
time being as unrestrictive as possible. Further discus- 
sion of the experimental foundations for th e assumptions 
can b e foun d in Zhytnikov & Nester (1994 ) an d in Wein- 
berg ( |1972| ) Misn er, Thorne & Wheefer~( ^973| ) and espe- 
cially WTTl (H993l §§2-3). 



For any such model, the important change in the formal- 
ism for the propagation of light through such a weak field 
metric is to change the Poisson equation relating the den- 
sity to the potential whose derivative is used to determine 
the bend angle of photons. The angular power spectrum 
of the convergence, n, can be written as an integral over 
the line-of-sight of the power spectrum of the density fluc- 
tuations (Kaiser 1992). For sources at a distance D s , 



£{l + l)C e /(2n) = f £ [n m H 2 D 2 ] 2 J ^t\l - tf 
s (k = t/D,a)~ 



A 



Nester 1994; Edery 199E; Kinney & Brisudova 2001; Uzan 



f 2 {k = l/D), (1) 



where t = D/D s , A^ aS8 (fc) = k 3 P(k)/(2TT 2 ) is the contri- 
bution to the mass variance per logarithmic interval physi- 
cal wavenumber and £(£+ \)Ci/(2-k) is the contribution to 
K 2 ms per logarithmic interval in angular wavenumber (or 
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equivalently multipole) I. The only change from the stan- 
dard result is that the Poisson equation relating the poten- 
tial to the density perturbations is modified from f(k) = 1 
to a functional form determined by the Poisson equation 
of the modified theory of gravity. On small physical scales 
(large wavenumber k), f(k) = 1 is required to be consis- 
tent with the known properties of gravity. 

If the sources have a range of redshifts then one simply 
integrates the above expression over the redshift distribu- 
tion of the sources. We shall assume throughout that 



^oc^cxphP/^) 4 ] 



(2) 



and fix 
uating Eq. 



the requirement that (z src ) = 1. In eval- 
, we will use the method of Peacock & 
Dodds ( |1996| ) to compute the non-linear power spectrum 
as a function of scale- factor. Throughout we shall u se the 
concordance cosmology of Ostriker & Steinhardt ( 1995| ) 
since it provides a reasonable fit to recent CMB, weak 
lensing and large-scale structure data. For this choice of 
parameters the lensing kernel peaks at z ~ 0.43 at a (co- 
moving) angular diameter distance of 1150/i _1 Mpc. 

In our calculation we only consider the propagation of 
rays through a known density distribution, and we model 
that known density distribution using a standard cosmo- 
logical model viewed as a means to interpolate the evolu- 
tion of structure with redshift. We do not attempt to self- 
consistently form the observed structures using the modi- 
fied gravitational potential]^. If we assume that all theories 
must match the local density distribution, the only conse- 
quence of this assumption is that the evolution of structure 

1 In the model described below, a linear fluctuation analysis sug- 
gests that long-wavelength modes would grow more slowly than the 
standard model would predict. Thus neglect of this effect is conser- 
vative if we start from an initially scale-invariant spectrum. 




Fig. 1. — The angular power spectrum, 1(1 + l)C^/(27r), vs. mul- 
tipole moment I for models with a = 1.0 and m = 0, 0.3, 1.0 and 
3. Oh Mpc" 1 . The sources are assumed to have (z src ) = 1- Spectra 
for other values of a can be roughly obtained by averaging the m = 
spectrum and the appropriate a = 1 spectrum (plotted here) with 
the relevant weights. 



implicit in Eq. (Q) uses the standard growth rates rather 
than those of the modified gravity. 

Examining the effects of modified gravity simply be- 
comes a question of considering different structures for the 
function f(k). In 4D, the metric, being symmetric, con- 
tains 10 functions. The 4 constraints of energy- momentum 
conservation reduce the number of free functions to 6. 
These 6 free functions can be decomposed under rotations 
as 2 scalar (density perturbations), 2 vector (vortical mo- 
tions) and 2 tensor (gravity wave) modes. Within the lin- 
earized theory there are a number of propagating modes, 
which have the form of Yukawa (exponential) potentials 



U(r; m)=G 



p(r')d 3 r' m \^f\ 



(3) 



Under a vari ety of reasonable assumptions Zhytnikov & 
Nester (1994) conclude that the most general metric de- 
scribes forces mediated by massive and massless s calar a nd 
tensor particles. We follow Zhytnikov & Nester (1994) in 
neglecting the vector modes, however we will allow arbi- 
trary couplings for the scalar and tensor modes. In general 
relativity in the weak field limit 



ffoo 

9ij 



(-1 
( 1 



2U) 
2U) % 



(4) 
(5) 



where U is the usual Newt onian potential. The metric of 
Zhytnikov & Nester ( 1994 ) has the same form, but with 
Yukawa potentials in addition to the Newtonian one. 

For test particles with i; « c or fluids with p -C pc 2 
only the time-time part of the metric is relevant, the con- 
tribution of the gij terms being suppressed by 0(v 2 /c 2 ). 
However, for light, the bend angle due to the potential is 
actually the arithmetic mean of the coefficients in goo and 
g^. Though the extra scalar and tensor modes can enter 
into the space-space and time-time part of the metric dif- 
ferently, we shall consider the 1 parameter family of models 
where these coefficients are equal. As Kinney & Brisu- 
dova ( 2001 ) discuss, the requirement that cluster mass 
estimates from galaxy dynamics, pressure equilibrium of 
the X-ray gas and gravitational lensing agree means that 
any modified gravity law must affect photon propagation 
in roughly the same was as it affects particle orbits. A 
modified gravity which differentially affects particles and 
photons will almost always lead to a discrepancy between 
these three cluster mass estimates. 

Thus in our model, in the weak field limit, the propaga- 
tion of light is the same as in standard general relativity, 
except that the potential is 



U(r) = (1 - a)U(f, 0) + aU{r, m) + 



(0) 



where • ■ • represents possible other terms of the same form 
as the second. We shall further simplify our calculation by 
considering only 1 correction term in what follows. In such 
a theory with one additional "field" , the function appear- 
ing in the estimate of the weak lensing power spectrum 
is 

/(fc) = (l-a) + a F ^ ?l (7) 

where a = for standard gravity and a ~ —0.9 and 
m ~ 50kpc in order to produce flat rotation curves 
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without dark matter (e.g. Sanders 1986). The correspond 
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Fig. 2. — The rms shear, smoothed with a 5' FWHM gaussian (top) 
or a 10' FWHM gaussian (bottom), predicted for the "concordance" 
cosmology with f2 ma t = 0.3, Ha = 0.7, h = 0.67 and erg = 0.9 as 
a function of a and m (h/Mpc). Contours are spaced every 0.001 
with bold contours indicating 0.005 (top), 0.01 and 0.015 (bottom). 
The stippled regions are consistent (at ltr) with the van Waerbeke 
et al. (2001) measurements. 



ing potential for an object of mass M simplifies to the New- 
tonian result, —GM/r, on small scales where mr <C 1, and 
has a different effective coupling constant, —GM(1 — a)/r 
on large scales, mr ;§> 1. 

Fig. |l| shows the anisotropy spectrum predicted for a 
range of models. If we limit the range of gravity (a > 0) 
then the shear fluctuations on large angular scales are sup- 
pressed, and if we extend the range they are enhanced. 
This should be a generic feature of any modification to 
the long-range force law. To obtain limits on the parame- 
ters in our model we calculated the rms shear expected in 
Gaussian windows with FWHM of 5' and 10' as a function 
of a and m (Fig. ^) . These predictions are consistent with 
the r ms she ar measured on these scales by van Waerbeke 
et al. ( 200l| ) only for models with parameters close to those 
of standard gravity. We can minimize the model depen- 
dence of the result by examining the ratio of the power 
at 5' and 10', as this largely removes any dependence of 
the result on the matter density and the normalization of 
the power spectrum. In Fig. we see that the data are 
consistent with standard gravity and a broad range of al- 
ternate theories. These theories are acceptable because 
our alternate gravity model has a 1/r potential on large 
scales so that when the 5' scale corresponds to a physical 
scale larger than to -1 , the change in the coupling constant 
a is degenerate with a change in the enclosed mass. For 
sources with a mean redshift of unity, the 5' scale corre- 
sponds to a length scale at the peak of the lensing kernel 
of approximately 1/t, -1 Mpc. 

Theories which do not return to a 1/r potent ial on large 
scales are relatively easy to rule out (see Walker 1994). As- 
suming that the bend angle of light remains proportional 
to the gradient of the projected gravitational potential, 
such theories predict that random lines of sight would be 
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Fig. 3. — The ratio of the rms shear on 5' and 10' scales for the 
same cosmology as the previous figure. Contours are spaced in steps 
of 0.05, increasing to top left. Thick contours are spaced every 0.5, 
starting at 1.5. The stipp led region is consistent (at la) with the 
van Waerbeke et al. (2001) measurements. 
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highly sheared and (de)magnifiedg in contradiction with 
observations. This problem can be traced to the lack of 
degeneracy between renormalizing the mass and adjusting 
the coupling consta nt. For example, ignoring the Kinney 
& Brisudova (2001) ansatz for permissible forms of alter- 
nate gravity, we could use the force law 



<j>'(r)/GM 



1 exp(— mr) 



rr 



(8) 



which is Keplerian for r < ro an d r » m" 1 but is a 1/r 
force law, producing a flat rotation curve, in between. The 
potential corresponding to this force law is 



cj)/GM = -l/r + Ei[-mr] 



(9) 



where Ei[x] is the exponential integral. The corresponding 
kernel for the weak lensing integral is 



/(*) 



1 



km — (k 2 + m 2 ) tan 1 (k/m) 
r a k(k 2 + m 2 ) 



(10) 



Figure ^ shows the angular power spectrum in this model 
for a range of scales ro and a large outer cutoff mT 1 = 
50/i _1 Mpc. Compared to normal gravity, the modified 
theories have enormously enhanced large scale power and 
very different shapes. 

3. DISCUSSION 

Current modified gravity theories tuned to explain the 
rotation curves of galaxies work in a standard cosmology 
because we measure rotation curves only where there are 

2 For example, for a logr potential and a Poisson distribution 
of lenses the convergence, k, of a source at D s (assumed to be 
much larger than the scale, ro, beyond which gravity is logr) is 



wag I 



^> 1 for any reasonable source density rt. 



10 



10- 




FlG. 4. — The angular power spectrum, i(£ + l)C(/(2n), vs. mul- 
tipole moment I for our second model with m _1 = 50h~ 1 Mpc and 
ro = oo, 20, 10 and bh' 1 Mpc. As ro — » oo the model becomes 
standard gravity. Note the change in shape and the enormous en- 
hancement in the power on large scales. 



baryons. We can see that the rotation curve is flat out 
to the limit where there are no more baryons to measure, 
but we cannot see that it is Keplerian as we approach the 
edge of the more extended dark matter distribution. If we 
could continue to trace rotation curves on larger scales we 
would see a growing difference between standard cosmo- 
logical models and theories using modified gravitational 
physics. 

Weak lensing allows us to do this experiment, although 
on such large scales we must sum over the contributions of 
all of the mass rather than consider the rotation curves of 
discreet objects. As we would expect qualitatively, increas- 
ing the strength of the gravitational field at long ranges 
predicts stronger weak lensing signals on large scales than 
standard cosmological models. Current measurements of 
the rms shear on scales of 5'-10' rule out the theories we 
consider in the parameter ranges where they could explain 
rotation curves without dark matter unless the deviation 
from normal gravity is limited to a restricted range of spa- 
tial scales from 10/i _1 kpc < r < l/i _1 Mpc. On larger 



scales the models must return to the 



force law of nor- 



mal gravity in order to be consistent with measurements. 

In standard cosmological models, once we postulate the 
existence of dark matter, the inner and outer scales appear 
naturally. On small scales the cooling of the baryons con- 
centrates the baryons relative to the dark matter and ren- 
ders them luminous and detectable. Thus, normal matter 
combined with normal gravity naturally explain dynamics 
on scales <C lO/i -1 kpc. On intermediate scales, dark mat- 
ter provides an additional source of density, which can be 
interpreted as an abnormal gravitational theory using only 
the visible baryons as sources. On large scales the universe 
returns to homogeneity, and the special properties of the 
1 jr 2 force law make the weak lensing power slowly dimin- 
ish on large scales. Abnormal, longer ranged theories lose 
the cancellation properties of the 1/r 2 force law on large 
scales, despite the increasing homogeneity of the density 
on these scales, leading to enormous enhancements in the 
strength of the weak lensing shear. Such strong shears 
are in gross disagreement with even the first generation of 
weak lensin g measurem ents o n these scales (van Waerbeke 
et al. |200C| ; Bacon et a l |200C| ; Kaiser et al |200C| ; Wittman 
et al. pOOOf ; Maoli et al. |2001t Rhodes et al. |2001fc van Waer- 
beke et al. [200 1| ). Thus, abnormal gravity theories must 
introduce two fine-tuning scales - an inner scale to explain 
flat rotation curves and an outer scale to force a return to 
Newtonian gravity on large scales - and these scales must 
coincidently match the scales produced by the dark mat- 
ter theory after evolving the universe for 10 billion years 
starting from initial conditions which are exquisitely de- 
termined from the cosmic microwave background. 

Finally, although we lack a formalism for estimating 
weak lensing in non- poten tial theories such as MOND, 
Mortlock & Turner (2001) have emphasized that weak 
lensing results should be generic, as it requires only that 
photons and particles have similar responses to gravita- 
tional fields. This similarity of behavior is observed on 
the relevant scales (Mpc) through the near equivalence 
of weak lensing, dynamical, and X-r ay de terminations of 
cluster masses (Kinney & Brisudova 2001). 
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